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Tensor network states are capable of describing many-body systems with complex quantum en- 
tanglement, including systems with non-trivial topological order. In this paper, we study methods 
to calculate the topological properties of a tensor network state from the tensors that form the 
state. Motivated by the concepts of gauge group and projective symmetry group in the slave- 
particle/pro jective construction, and by the low-dimensional gauge-like symmetries of some exactly 
solvable Hamiltonians, we study the d-dimensional gauge structure and the d-dimensional symmetry 
structure of a tensor network state, where d < d spa ce with d spa cc the dimension of space. The d- 
dimensional gauge structure and d-dimensional symmetry structure allow us to calculate the string 
operators and d-brane operators of the tensor network state. This in turn allows us to calculate many 
topological properties of the tensor network state, such as ground state degeneracy and quasiparticle 
statistics. 



I. INTRODUCTION 

Topological order [1-3] is a new kind of ordering in 
many-body quantum states. It represents new states of 
quantum matter beyond the symmetry breaking states 
characterized by Landau symmetry breaking theory. 
This new kind of order and the new states of matter as- 
sociated with it open up a whole new research direction 
in condensed matter physics. [4-18] Intuitively, topologi- 
cal order corresponds to a pattern of long range quantum 
entanglement in the ground state. 

Traditional mean-field approaches to calculating the 
T = phase diagram of a quantum system are based 
on direct product states which have no long range quan- 
tum entanglement. These approaches seem to exclude 
topological states from the very beginning. To study a 
phase diagram that may contain topological states, we 
have to find a way write down states with non-trivial 
long range entanglement. The tensor network state 
(TNS) [19-22] method is one way to produce states with 
long range entanglement. [23-25] We can develop a mean- 
field approach for topologically ordered states based on 
TNS. [26-31] 

Physically, topological order is characterized and de- 
fined through measurable quantities, such as ground state 
degeneracy on a torus or other compact space [1, 32] and 
fractionalized quantum numbers and statistics of quasi- 
particles [33, 34]. If we claim that TNS can describe 
topological states, then it is natural to ask how to calcu- 
late those measurable topological quantum numbers from 
the TNS representation. In this paper, we will address 
this question. We will show how to calculate measurable 
topological properties of a TNS. 

Our approach is motivated by the gauge group and pro- 
jective symmetry group in the slave-particle/projective 
construction [35, 36] and by the low-dimensional gauge- 
like symmetries of some exactly soluble Hamiltonians. [37] 
We will introduce the local gauge and local symmetry 
structures of a TNS from the tensors in the TNS. This 
allows us to obtain the string operators (or more gener- 
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FIG. 1: Tensor-network - a graphical representation of the 
tensor network wave function (1), (a) on a ID chain or (b) 
on a 2D square lattice. The indices on the links are summed 
over. 

ally, the d-brane operators) of the TNS, which in turn 
allows us to calculate the topological properties of the 
TNS. 



II. TNS, IDEAL WAVE FUNCTIONS, AND 
IDEAL HAMILTONIANS 

What is a TNS? As an example, let us consider a TNS 
on a square lattice, where the physical states living on 
each vertex i are labeled by mi. The TNS is defined by 
the following expression for the many-body wave function 

*t(W) = E ^ejifT™? hkg TZU T Zhi ■ ■ (!) 

ijkl--- 

Here T™^^ is a complex tensor on vertex i with one 
physical index rrii and four inner indices i,j,k,l,---. The 
physical index runs over the number of physical states 
Dp on each vertex, and the internal indices runs over D 
values. Note that tensors on different vertices can be 
different. The TNS can be represented graphically as in 
Fig. 1. 

We see that a TNS is characterized by tensors T^-^ 
and by a network specifying the connectivity. Can we 



FIG. 2: (Color online) The patch Pi and the patch Pj. Pi 
and Pj only differ by a displacement. 

calculate the topological properties of the many body 
state ^rdiTii}) from its defining tensors T^ kl 7 To an- 
swer this question, we need to define the problem more 
completely by introducing a Hamiltonian Ht such that 
^T{{m,i}) is the exact ground state of Ht- We will call 
Ht an "ideal" Hamiltonian and ^T({mi}) an "ideal" 
wave function. 

The Hamiltonian Ht has the following local form 

H T = J20 Pi (2) 

i 

where Op t is an operator that acts only on states in the 
patch Pi (see Fig. 2). How can we choose the operator 
Op i so that v E , T({wj}) is, hopefully, the only ground state 
of iJp? In the following, we will describe a construction 
proposed in Ref. 38. 

Let pp i be the density matrix of the state ^T({mi}) 
on patch Pi obtained by tracing out the physical degrees 
of freedom outside Pi : 

p Pi =Trp.|* T )(* T | (3) 

where Pi represent the region outside of Pj and Tr^i is the 
trace over all the states in region A. From the structure 
of the TNS, we can deduce that pp i must have the form 

p Pi = E MfeXM ( 4 ) 

o=l 

where \tjj a ) and \4> a ) are states on the patch Pi, and Lp i is 
the number of links on the boundary of the patch Pi (for 
example Lp i = 10 for the patch in Fig. 2). Equivalcntly, 
pp i must have the form 

N P . 

PPi=52pa\<Pa)(Va\,- N Pi <D 2L ^, (5) 
a=l 

where \ip a ) are normalized states on the patch Pi. This 
means that pp i can be constructed from only Np i < 
D 2Lp i states on patch Pi. Since the number of states 

aL 2 

on Pi grows like D p Pi , a = 0(1) for a large patch (here 
rfspace = 2), the following projection operator 

N P . 
a=l 



FIG. 3: A graphical representation of a gauge transformation. 




FIG. 4: (Color online) A 1-brane and a 2-brane. 
is non-trivial in the large patch limit {i.e. 1 — Pp ± ^ 0) 

aL 2 

because N Pi < D 2Lp i < D p Pi . If we choose 

H T = ^ l - p P^ ( 7 ) 

i 

the TNS ^t will be the exact ground state of Ht- We 
hope that in the large patch limit, the projection Pp i is 
highly restrictive and 'Fp is more or less the only ground 
state. 

Now the problem can be better defined. From the ten- 
sors TY\- kl , we can determine an "ideal" Hamiltonian Ht 
(7). Assuming H T is gapped, the topological properties 
of the TNS are just those of the topological phase of 
Ht- However, in this paper we will not try to calculate 
the topological properties by directly diagonalizing Ht- 
This task is generically difficult because [Pp i ,Pp j ] ^ 0. 
Instead, we will find a way to calculate the topological 
properties from the structure of the tensors T™^ kl di- 
rectly. 

III. LOCAL GAUGE STRUCTURE OF A TNS 

One of the important properties of a TNS is that dif- 
ferent tensors TJpLjy can give rise to the same physical 
wave function $t({™i})- To be more precise, two sets 
of tensors T™ i - kl and T^- kl related by a "gauge transfor- 
mation" (GT) 

A-iTiJ'r'd'u' (8) 

= ^ t (Ai — X ,i)l' l(Ai+ x f i) r i r (Ai+ y _i) u i u {Ai—y t i) d> dTi^rdui 
Irud 

Yj^y^.'-\,f.^ a ;j ,\, (9) 

fc 
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give rise to the same wave function ^T({mi\) up to some 
scaling factors (see Fig. 3). In the above Ai and Ajj are 
scaling factors that do not depend on the indices i, j,k, .... 
Here Aij and Aj t i are D x D invertible matrices de- 
fined on the links A GT is called a "d-dimensional 
gauge transformation" (d-GT) if Aj.i ^ 1 only on links 
within with a "d-brane" . By "d-brane" we mean a region 
near a d-dimensional subspace (i.e. a d-dimensional sub- 
space with a finite "thickness" £ as illustrated in Fig. 4). 
Note that the "thickness" £ is fixed for a TNS, regardless 
the size of the d-dimensional subspace and the value of 
d. The size L of the d-dimensional subspace is typically 
much bigger than £: $/L4 0. 

Using d-GT we can introduce the concept of a "d- 
dimensional full invariant gauge group" (d-fIGG) as in 
Ref. 35. d-fIGG is a property of the tensors T^\- kl and a 
fixed d-brane. d-fIGG is a group formed by all the d-GT 
in the d-brane that leave T^ n i - kl invariant up to a scaling 
factor. 

We note that if d-fIGG is non trivial then we can use d- 
GT in d-fIGG with different d-branes to construct a non- 
trivial d'-flGG for all d' > d. Thus, to reveal the new 
gauge structure that appears at every dimension d, we 
will introduce the "d-dimensional invariant gauge group" 
(d-IGG). The d-fIGG contains a normal subgroup whose 
elements are in (d — l)-fIGG. The d-IGG is defined as 

d-IGG ee d-fIGG/(d - l)-fIGG (10) 

The d-IGG for the tensors T t ™ . feJ defined this way may 
depend on the topology of the d-brane. 

As an example of this notation, we note that for a 
dspaco-dimensional TNS, the d spa ce-IGG is like the group 
of uniform global gauge transformations in a gauge the- 
ory. 

IV. LOCAL SYMMETRY STRUCTURE OF A 
TNS 

In the previous section we discussed gauge transforma- 
tions that leave the tensors of the TNS invariant (up to a 
scaling factor). These gauge transformations define the 
d-fIGG. In this section, we will discuss gauge transfor- 
mations that do change the tensors TV n, i - kl . 

Let us first introduce some concepts. A "local physical 
transformation" of a TNS is generated by unitary D p x 
Dp matrices Si, TOm ' acting on the physical indices m (here 
Dp is the range of the physical index m): 

T™ijkl ~* ^ Si,mm'Ti?ijkl (H) 
m' 

A "d-dimensional local physical transformation" (d-lPT) 
is a "local physical transformation" that is non-trivial 
(i.e. Si_ mm i =/= 5 m m') only on a d-branc. A "d- 
dimensional symmetry transformation" (d-ST) of a TNS 
is a d-lPT that leave the TNS invariant. All the d-ST 
on a fixed d-brane from a group which will be called the 
"d-dimensional full symmetry group" d-fSG. 



The d-fSG is very similar to the "low-dimensional 
gauge-like symmetries" introduced in Ref. 37. A differ- 
ence between the two concepts is that d-fSG is a symme- 
try of a ground state while the "low-dimensional gauge- 
like symmetries" are symmetries of an exactly solvable 
Hamiltonian. As pointed in Ref. 37, the low-dimensional 
gauge-like symmetries are useful tools to calculate the 
topological properties of a model Hamiltonian. Simi- 
larly, d-fSG should be a useful tool to calculate topo- 
logical properties from a ground state. 

We note that if d-fSG is non trivial, then we can use 
elements in d-fSG on different d-branes to construct a 
non-trivial d'-fSG for all d' > d. Thus, to reveal the new 
structure that appears in each dimension d, we introduce 
the "d-dimensional symmetry group" d-SG. As in the 
previous section, the d-fSG contains a normal subgroup 
that is formed by elements in (d— l)-fSG. We define d-SG 
as 

d-SG ee d-fSG/(d - l)-fSG (12) 

We can make contact with the string-net framework 
by noticing that the elements of 1-SG are actually the 
ideal string operators [12, 16, 39] that satisfy a "zero 
law". [40] Similarly, the elements of d-SG are the ideal 
d-brane operators that satisfy a "zero law". [41, 42] So 
d-SG and "low-dimensional gauge-like symmetries" are 
closely related to the string operators and the d-brane 
operators that play a key role in topologically ordered 
states. [14-16] 

The above discussion of d-SG (or "low-dimensional 
gauge- like symmetries") is for a generic topological state. 
Now the question is how to calculate d-SG for a TNS? 
Naively, we would like to say that d-SG are formed by lo- 
cal physical transformations that leave the tensors T^\- kl 
invariant. But due to the gauge structure of the TNS, 
the invariance of the TNS and the invariance of the ten- 
sors TV"-- kl are not equivalent. Instead, d-SG are formed 
by local physical transformations that leave the tensors 
T[}jki invariant up to gauge transformations. Thus the 
gauge transformations that change the tensors play an 
important role in calculating d-SG for a TNS. 

In the above, we have defined a notion of d-IGG and 
d-SG for a TNS. We will sec that both d-IGG and d- 
SG play a very important role in characterizing a TNS. 
Many (if not all) topological properties of a TNS can be 
calculated from its d-IGG and d-SG. In the following, we 
will carry out the this program for the simplest TNS with 
a non-trivial topological order. 



V. Z 2 TOPOLOGICALLY ORDERED STATE 

A. Z2 String Condensed State and its Ideal 
Hamiltonian 

The simplest topological phase that we can consider 
is the Z 2 topological phase. [8, 9] Such a phase can 
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These operators can be taken to be 



FIG. 5: (Color on line) A Kagome lattice viewed as the links 
of a honeycomb lattice. 
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FIG. 6: (Color online) A spin configuration of <r z = 1 (filled 
dot) and a z = — 1 (open dot) on the Kagome lattice can be 
viewed as a string configuration on the honeycomb lattice. 
The end of an open string costs an energy U. The closed- 
string sector has an energy gap g. The dashed lines describe 
the dual lattice of the honeycomb lattice. The thick blue lines 
are strings in the honeycomb lattice. The thick green lines are 
strings in the dual lattice. 



be realized through a spin- 1/2 model on the Kagome 
lattice: [12, 16] 

= 17^(1 -Q z ) -Bp) 
/ P 

Qi= n ^ b p= n <■ 

legs of J edges of p 

Here we have viewed the Kagome lattice as the links of 
the honeycomb lattice (see Fig. 5). The vertices of the 
honeycomb lattice are labeled by J, the links (which are 
the sites of the Kagome lattice) by i, and the faces by 
p. J^i sums over all vertices and J2 P over an faces. The 
Hamiltonian (13) indeed has the form of a summation of 
projectors as in eqn. (7). 

We can interpret a given configuration of spins in terms 
of strings by viewing the state a z = 1 as the absence of a 
string and the state a z = — 1 as the presence of a string 
(note the choice of a z here). The above Hamiltonian has 
the property that, in the large U > limit, the low energy 
Hilbert space is the space of all spin configurations that 
contain only closed strings (see Figure 6). For g > 0, the 
ground state of H is the equal weighted superposition of 
all closed string states |0) = ^xciosed 1^0- ^ * s useiu l to 
consider string operators that add a string to the system. 



x{c) = n of 



(14) 



where C is a curve running along the edges of the lattice. 
The operator X(C) is a string creation operator because 
a x acts as a spin flip operator for eigenstates of a z = m 
so that we are simply flipping the spins (adding a string) 
along the curve C. Note that the operator X(C) can also 
remove a string because the string is its own anti-string 
meaning that two strings can annihilate each other. 

The ground state |f2) has the remarkable feature that it 
is an eigenstate of X(C) for all C with eigenvalue 1. We 
say that the strings created by X(C) have condensed in 
the state Another useful point of view comes from 
thinking about the Z2 condensed state in terms of Z2 
gauge theory where the operator X{C) is nothing but a 
Wilson- Wegner loop operator. In the string condensed or 
deconfined phase the Wilson- Wegner loop X{C) satisfies 
a "zero law" as opposed to an area or perimeter law. [40] 
This "zero law" expresses the existence of strings at all 
scales in the ground state. 

For an open string operator that satisfies the "zero 
law" , the action of the string operator will create a pair 
of excitations at its two ends. Such excitations will in 
general have fractional statistics and fractional quantum 
numbers. For the X(C) string, the excitations at its ends 
correspond to Z2 charge. Thus we will call the X(C) 
string an electric string. 

There is another string operator that we can consider, 
which we define as 



Z(C*) = J] at 



(15) 



where C* is a curve along the links of the dual lattice (see 
Figure 6). It is characteristic of the Z2 string condensed 
phase that these strings are also condensed: Z(C*)\Q) = 
|0) for any loop C*. 

For the Z(C*) string, the excitations at its ends corre- 
spond to Z2 vortices. Thus we will call the Z(C*) string 
a magnetic string. 

This was all for the case of a system on an infinite 
plane where there is a unique ground state. Things 
change qualitatively when the system lives on a manifold 
with non-trivial topology. When the manifold in question 
has non-trivial 1-cycles (non-contractible loops) then the 
string operators corresponding to these non-contractible 
loops become interesting observables. We will return to 
the case of non-trivial topology later. 



B. Z2 Tensor Network Representation 

The Z2 condensed state, and many other topologi- 
cally ordered states, can be systematically represented in 
terms of a TNS.[24, 25] Fig. 7 and 8 represent the tensor 
network on the honeycomb lattice graphically. Note that 
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FIG. 7: A tensor network formed by vertices and links. The 
links that connect the dots carry indices a,b, .... Each triva- 
lent vertex represents a T-tensor. The vertices on the A- 
sublattice (open dots labeled by I a) represent T abc . The ver- 
tices on the B-sublattice (shaded dots Ib) represent T' abc . The 
dots on the links represent the <? m -tensor g™ b . In the tensor 
trace tTr® i A T®i B T' (3ig mi , the a, b, c, ... indices on the links 
that connect the dots are summed over. The tensor trace pro- 
duces a complex number $(mi,m2, ...) that depends on m ; 
which can be viewed as a wave function. 




a 



A » J - 

m 
(b) 



FIG. 8: The graphic representation of (a) the T-tensor, T abc , 
and (b) the g m -tensor gT A b B - 



the tensor network is formed by three tensors T abc , T' abc , 
and g™ b . To describe the Z2 condensed state, we choose 
D = 2 and 



ffoo 

Tabc 
T a bc 



= 1, 

rpl 

abc 
rpl 
- 1 - abc 



9u =1, other 3 ^ = 0; 
= 1 if a + b + c = even, 
= if a + fo + c = odd. (16) 



where a, b, c = 0, 1. 

In the following, we will calculate the string operators 
directly from the above tensors using the d-IGG and d- 
SG of the tensor network. This allows us to understand 
the topological properties of the TNS directly from its 
defining tensors. 



C. d-IGG of the Z 2 Tensor Network 

What is the d-IGG of the above Z2 tensor network? 
First, we find that the 0-IGG is formed by the following 



gauge transformations: 

(An)ab = Ai/(5 ab , (Au) ab = XnSab- (17) 

These transformations only change the tensors g, T and 
T' by a scaling factor. 

Next, let us calculate 2-IGG. We consider the "gauge 
transformations" that leave g"l on link-i invariant, up to 
scaling factor: 

K ig % b , =Y,{An) a >a{Aji) blb g™ b , form = ±l. (18) 

ab 

Note that J and J label the vertices of the honeycomb 
lattice and i is the link that connect the two vertices I 
and J. The above equation requires that An) and Aji) 
be diagonal. Using the "local gauge transformation" : 



9ab -> ^i^'a A 2 <V6SS, 

ab 

we can fix (An) o = (^ji)oo = 1 and 
{A H )u(Aji) n = 1. 



(19) 



(20) 



From cqn. (9), we find that An must be also diagonal. 
Using the 0-GT in 0-IGG=0-fIGG we can transform the 
diagonal An into the following form 

(Aa)oo = 1, (Auju = l/(A H )u. (21) 
An must leave T invariant and thus 

Ta'b'c' = ^ (Aii) a / a (Aji)b'b(Aki)c'cT a bc, (22) 

abc 

where i, j, and k are the three links that connect to the 
vertex I. The above gives us 

(^1)11(^1)11 = (Ajj)ii(Afei)ii = {AkiJuiAuJu = 1. 

(23) 

Such a non-linear equation gives two sets of solutions 

A a =An = l, (24) 

and 

A iI = A Ii =a z . (25) 

where we have used eqn. (21). So 2-IGG of the Zi tensor 
network contains only two elements given by the above 
two expressions. Such a 2-IGG is a Zi group. 

Based on this discussion, we see that we can ob- 
tain the 2-IGG from the 2-fIGG by removing the 0-GT, 
i.e. 2-IGG = 2-fIGG/O-fIGG. This implies that the 1- 
IGG is the trivial group containing only the identity. 



6 




FIG. 9: (Color on line) The Z — 2 gauge transformation An — 
An = a z acts only on the shaded region. The tensor g is 
changed to g on the boundary of the stripe. 

D. Ground State Degeneracy from d-IGG 

Using the ci-IGG of the TNS, we can calculate some 
topological properties of the TNS. In the Z 2 condensed 
state we found that 2-IGG=Z2. Now we will show that 
the ground state degeneracy of the Z 2 condensed state 
on a torus is the same as the ground state degeneracy of 
Z 2 gauge theory on a torus. 

If we apply the Z 2 gauge transformation An = An = 
<j z every where, the tensors T, T" and g are not changed. 
If we apply the Z 2 gauge transformation only on a stripe 
(see Fig. 9), then the tensors on the boundary of the 
stripe will be changed. For the particular choice of the 
stripe in Fig. 9 the tensor g is changed to g on the bound- 
ary of the stripe, where 

3o"u =1 = l, SIT' 1 —!, other 5™ b = 0. (26) 

We note that the change g —> g along the boundary of 
the stripe is generated by the Z(C*) string operator (see 
eqn. (15)). Therefore the Z 2 gauge transformation on a 
finite region defines a type of string operator along the 
boundary. 

Since the application of the Z 2 gauge transformation 
does not change the energy of the state, we find that the 
application of the Z(C*) string operator also does not 
change the energy of the state. We note that the action 
of two Z(C*) string operators along the top and the bot- 
tom boundary of the stripe in Fig. 9 is equivalent to a Z 2 
gauge transformation on the stripe. In other words, the 
tensors are invariant under the action of the two Z(C*) 
string operators, up to a Z 2 gauge transformation on the 
stripe. Thus the product of the two Z{C*) string opera- 
tors is an element in 1-SG. The action of the two Z(C*) 
string operators does not change the state. 

On the other hand the action of one Z(C*) string op- 
erator is not equivalent to a Z 2 gauge transformation and 



is not an elements of 1-SG. Such an action will change 
the state, but not the energy. Applying such a Z(C*) 
string operator to the Z 2 condensed state on torus along 
the x- and y-directions will give us the four degenerate 
ground state of the Z 2 condensed state on a torus. This 
is how do we calculate the ground state degeneracy from 
the tensors in the TNS. 

We would like to stress that the Z(C*) string opera- 
tor is directly obtained from the Z 2 gauge transforma- 
tion. The Z 2 group structure of the gauge group also 
determines the Z 2 structure of the Z(C*) string opera- 
tor: [z(c*)\ 2 = i. 



E. Z2 d-SG and the Electric String Operator 

We have seen that the d-IGG of a TNS can be calcu- 
lated by finding all the gauge transformations that leave 
the tensors invariant (up to a scaling factor). Similarly, 
the e?-SG of a TNS can be calculated by finding all the 
combined local physical transformations and the gauge 
transformations that leave the tensors invariant. 

For the case of the Z 2 condensed state, we have 
seen that the magnetic string operator Z(C*) on a con- 
tracture loop C* is a 1-ST (an element of 1-SG). The 
tensors in the tensor network arc invariant under Z(C*) 
followed by a 2-GT on the disk D bounded by C* . The 
1-SG generated by such a 1-ST, Z(C*), will be called a 
non-local 1-SG. 

There is another type of 1-ST, such that the tensors in 
the tensor network are invariant under the 1-ST followed 
by a 1-GT on the same curve on which the 1-ST is de- 
fined. The 1-SG generated by this type of 1-ST will be 
called a local 1-SG. 

The Z 2 condensed state also has a local 1-SG, which 
is generated by the electric string operator X(C) (see 
eqn. (14)). One can directly check that the tensors g, 
T, and T" on the tensor network are invariant under the 
action of X(C) followed by a gauge transforation on the 
loop C with 

An = A n = a x (27) 

on the loop. 

We have seen that the ground state degeneracy on a 
torus can be calculated from the string operators. We 
can also calculate the quasiparticle statistics from the 
string operators. This is because an open string operator 
creates a pair of quasiparticles at its ends. The electric 
string operator X(C) creates a pair of Z 2 charges, and 
the magnetic string operator Z(C*) creates a pair of Z 2 
vortices. The open string operators can also be viewed 
as hopping operators for the quasiparticles. The statis- 
tics of the quasiparticles can then be calculated from the 
algebra of these hopping operators {i.e. the open string 
operators). [14] 
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VI. DISCUSSION AND CONCLUSIONS 

In this paper, we viewed a TNS as an ideal wave func- 
tion that is the exact ground state of a ideal Hamilto- 
nian. When the ideal Hamiltonian has a finite energy 
gap, the TNS can represent a quantum phase with non- 
trivial topological order. We argued that the topological 
properties of such a phase can be calculated just from 
the tensors that form the TNS. 

Motivated by the gauge group and projective symme- 
try group in the slave-particle/projective construction 
[35] and by low-dimensional gauge-like symmetries of 
model Hamiltonians[37], we introduced the d-IGG and 
the d-SG for a TNS. Using the d-IGG and the d-SG of 
a TNS, we can calculate the string operators (or d-branc 
operators) of the TNS, as we demonstrated in the simple 
Z 2 condensed state. Many topological properties of the 



TNS, such as ground state degeneracy and quasiparticle 
statistics can be calculated from the resulting string op- 
erators. This allows us to identify the topological order 
of the TNS just from the tensors that form the TNS. 

While completing this paper, we became aware 
a preprint by Schuch, Cirac, and Perez-Garcia 
(arXiv: 1001. 3807) where a class of TNS satisfying the "G- 
injective" condition are studied (injectivity means that 
one can achieve any action on the inner indices a, b, c, ... 
by acting on the physical indices m). From the symme- 
try properties of the tensors (similar to our d spaco -IGG), 
they can also calculate ground state degeneracy and other 
physical properties from the tensors of the tensor net- 
work. 

This research is supported by NSF Grant No. DMR- 
0706078 
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